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We get the distributions of shower partons initiated by heavy quarks c and b by 
studying the fragmentation functions in the framework of the recombination model. 
The transverse momentum spectra of heavy flavored mesons are predicted with these 
distributions. We find that the contribution from the recombination of thermal- 
shower partons is an important part in the total spectrum for the mesons. We predict 
the heavy flavored meson productions for different centralities with the heavy quark 
fugacities fitted by the experimental data of J/ip transverse momentum spectra in 
Au+Au collisions. 
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I. INTRODUCTION 



Heavy quark production in relativistic heavy ion collisions is a subject of current interest 
for understanding the quark-gluon interactions. Since heavy quarks are produced in the 
initial hard scattering processes at RHIC, they are sensitive to probe the medium formed 
in the collisions as they may lose energy by gluon radiation during propagating through 
the medium. Much progress has been made to improve theoretical understanding of the 
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production mechanism of heavy mesons containing a heavy quark (c or b quark), such as the 

.,„„,....„..,.n...,,... .„.„,.nT, 

Heavy Quark Effective Theory (HQET) [2| based on QCD. But the techniques of HQET 
can not be apphed to Bc(bc) meson since the charm quark is not sufficiently hght. There 



are two mechanisms investigated for the production of 6-flavored mesons: recombination 
and fragmentation. And it has been pointed out 14 1 that the dominant mechanism is due to 

n 

fragmentation, especially in the region of large transverse momentum. In Ref.[5[, we have 
shown the important contribution from the recombination of thermal-shower partons to the 
transverse momentum spectrum of J/ip in central Au+Au collisions. 

Up to now, the recombination model for particle production has achieved a great success. 
Especially the fragmentation process has been explained as a recombination of semi-hard 
shower partons in a jet. The model displays how the final state particles are produced 
from the recombination of thermal and shower partons created by hard partons. The de- 
termination of the shower parton distributions (SPD) has been obtained by studying the 
fragmentation functions (FF) in the framework of the recombination model p, |6| . With the 
parametrized SPDs, the model has offered good agreement with the experimental data on 
the hadron spectra from low to intermediate |7|, and successfully explained the Cronin 
effect on pion production in d+Au collisions js]. Also, based on the good description of light 
hadrons productions, the model has been extended to the productions of strange hadrons 9| 
and J/ip 5^ in central Au+Au collisions. In this paper we complete the SPDs initiated by 
hard partons of c and b quarks and furthermore predict the transverse momentum spectra 
for heavy-fiavored mesons. 

This work is organized as follows. We first review the basic formulism for the meson 
production in the recombination model in the next section. Then we show the results of 
charmed mesons and 6-fiavored mesons productions in central Au+Au collisions at ^/sNN = 
200 GeV in Section 3 and 4, respectively. In Section 5 we predict the heavy fiavored meson 
transverse momentum spectra for non-central centralities. A brief summary is given in the 
final section. 
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II. MESON PRODUCTION IN THE RECOMBINATION MODEL 

The contribution to the transverse momentum spectra of meson consists three terms 
because there are two components of parton sources: thermal (T) partons and shower (S) 
partons which are originated from hard partons. In this paper we substitute p for pt to 
denote the transverse momentum of the meson. Then the production can be expressed as 
the sum of three terms 

dNM d{Nl/ + Nlf + Nff) ^ 



(Pp d?p 



Ref.[lO|] provides a theoretical description of the recombination process in a dense parton 
phase, where the pure thermal contribution (TT) to meson spectrum can be written as 

^ = C^mM^^27.7./o[^^] [ cix|0M(x)pfcM(x,p). (2) 

Here 

kM{x,p) = M^^^{^ml+pl + ^ml + pl% (3) 

with the momenta of the two constituent quarks pi = xp and p2 = (1 — x)p. Iq and Ki 
are the modified Bessel functions. Mt is the transverse mass of the meson and At = PqIT 



is the transverse area of the parton system with the radius po = 9 fm 10|. It was assumed 



that hadronization occurs at r = 5 fm with temperature T = 175 MeV in the parton phase 

n 

[10|, which is consistent with predictions of the phase transition temperature at vanishing 
baryon chemical potential from lattice QCD Uj . The transverse flow rapidity t]t is defined 



by a flow velocity with vt = tanhr?T- In and ji, stand for the fugacities of the constituent 



quarks a and b, respectively. In Ref. [10| the fugacities of light quarks are 7u = 7d = 1 and 
for strange quarks 7^ = 0.8. Every quark has 3 color and 2 spin degrees of freedom. So we 
use the meson degeneracy factor Cm = (3 x 2)^. 

(pAiix) in Eq.(l2]) is the wave function of the meson in the momentum space. We can 
describe the wave function from the deflnition of the recombination function (RF). The RF 
for the process of qq' — > M with quark q at momentum fraction Xi and q' at X2 to form a 
meson at x is expressed as [12|] 

rm{x„x2,x) = -J— fi^y r^^v^i^ + ^ - 1), (4) 



B{a, b) \ X J \ X J X X 
where B{a, b) is the beta function. The RF Rm is an invariant distribution which is related 
to the non-invariant probability 6*1^(2/1,1/2) of flnding the two constituent quarks in M with 



momentum fractions yi and y2 through [jj 

RMixi,x2,x) = yiy2GM{yi,y2),yi = Xi/x. (5) 

The normahzation factor 1/B{a,b) in Eq.(jll) is determined by the following requirement 

/ dyi dy2GM{yi,y2) = I- (6) 

Jo Jo 

In fact, Guiyi^yi) relates with |0m(2/i)P by Ga/ ?/2) = |0Af(z/i)p5(2/i + 2/2 - !)• Conse- 
quently, the wave function of a meson in the momentum space corresponding to Cm (1/171/2) 
can be written as 

\'^M{yi)? = ^^^yt\^-yir'- (7) 

For a meson, the ratio of the average momentum fraction is proportional to the masses of the 
constituent quarks a/h = Xi/x^ = ma/mb^^, i.e., for meson, rric — 1.5 GeV, ruu d ~ 0.3 
GeV, a/b = rriu/mc — 1/5. 

The first step of calculating the thermal-shower recombination term {TS contribution) is 
to get the shower parton distribution. Using the parametrized distribution S{{z) [5, 6|, we 
can determine the distribution of shower parton ?' with transverse momentum p in central 

n 

Au+Au collisions as [16[ 

dq 
1 



where 



S{p) = T.j jF,{q)Si{p/q), (8) 



with //(fc) = fi{k) ■ {2Tif/E [Sj. The distribution fi{k) = dNj^^'-'^/d^kdy of hard parton 
, ju=t afte. ha.d Au+Au co.U.ons at ^ ^ 200 GeV at „nd-.apid.ty can 

be found in Refs.[l^ and [15|]. /3L is the explicit dynamical medium factor to describe the 
energy loss effect with PL = 2.39 for i =light quarks, gluon in Au+Au collisions for — 20% 



centrality 16| . Recent measurements of the transverse momentum distributions and nuclear 



modification factors of non-photonic electrons from heavy quark decays at high px show a 



suppression level of heavy quarks similar to light quarks 
loss factor 6L for i = c,b as that for light quarks. 



17|. So we choose the same energy 



In Ref. 



10|], the TS term can be calculated as 



d'p "Vs ^ (27r) 

Ua{R,Pl)Sb{R,P2) + Sa{R,P^)uJk{R,P2)l (10) 
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with the phase space distribution u{R,p) of thermal partons. u^{R) is the future oriented 
unit vector orthogonal to the hypersurface S defined by the hadronization volume. Because 
the thermal parton distribution density with momentum p is defined as jl^ 

91 



Tip) 



(27r)3 ^(-^'P)' 



(11) 



with g = Q coming from the color and spin degeneracy of a quark, we can rewrite Eq. (JTOll as 



dNlf 



(Pp 



Cm / dx\(t)M{x)\' 



.Ta{pi)Sh{p2) , Sa{pi)Tb{p: 



+ 



(12) 



91 aX 91b{^- x) 

We get the thermal parton transverse momentum distribution density from the thermal 
parton spectrum given in Ref. 10| by 



Tip) 



dp^dy 



T 



(13) 



The shower-shower recombination term {SS contribution) is equivalent to the FF 16 1 

^9 ^^,t.Pr^MtP^ 



pdp 



f '-^FM'-Di 

V ~^ J Q Q 



pOp 



where 



(3L 



0L 



dkkfi{k), 



(14) 



(15) 



and D^^ is the FF of quark i splitting into meson M. 

In the following sections, we show the results of heavy flavored meson transverse momen- 
tum spectra. 



III. PRODUCTIONS OF CHARMED MESONS 

Anomalous J/ip suppression is regarded as an important tool to probe the hot dense 
matter or quark-gluon plasma (QGP) created in heavy ion collisions at RHIC energies 18 |. 
The study of open charmed mesons {D and D) may help to explain the suppression of 
charmonium states {J/ip, Xc, In Ref. jsl, we have obtained the fugacity of c quark 

7c = 0.26 and Vt = 0.3c by fitting the measured J/ip transverse momentum spectrum in 
the region of low p^. Recently STAR collaboration have measured production at low 
transverse momentum (px < 2 GeV/c) in Au+Au collisions at RHIC energy [l9|]. To fit the 
experimental data with the thermal-thermal term of Eq.([2]), the factor of vt is adjusted to 



be Vt = 0.42c, which is a little larger than the fitted results in Ref.|5[, where vt = 0.3c, 
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Since Jjip is produced at the early stage in the heavy ion colhsions, and is produced 
later, it's reasonable that the radial flow velocity vt of is lager than that for J /ip. 

In Ref. jsl, the SPDs initiated by hard partons of charm quark and gluon have been 
achieved at flxed starting scale Q = 2mc = 3 GeV/c. Other SPDs with i = q,q,g and 
j = q, s, q,s{q = u, d, s) denoted by K^s, L, G, Lg and Gg in Ref. (6|] are obtained at Q = 10 
GeV/c. In order to set the SPDs above at the same scale, we investigate the SPDs K^s etc. 
at Q = 3 GeV/c with the method provided in Ref. 20|]. In calculating the thermal-shower 
{TS) term with Eq. (fT2|) . there are two kinds of shower partons. The sources of shower 
parton distributions of c and u given by Eq.(l8]) are different. For charm shower partons the 
dominated contributions are from i = c,g 5|], while for u shower partons they can come from 
^ = q^q^g- other contributions to the involved shower partons are neglectable because of 
either lower contents in the shower or less number of initial hard partons. So the summation 
formula Eq.(|8]) of u shower parton distribution in Au+Au collisions can be expressed as 



^ip) = I ^ii^Nsi^J + Li-))FM + Li^){FM + Uq) + F-M)) + (16) 
J q q q q q 



For D^{cu) meson, the RF is determined by Eq.(|l]) with a/h = l/5(a = 1) as mentioned in 
Section 2. 

In Eq.dH]) i = g,c because the D° FFs of u,d,s quarks are assumed to be zero at the 



starting scale 2l|. With the known FFs of D^" j5| and D^" [2l], we can calculate the 
shower-shower [SS) term. The results of production are shown in Fig{T] In order to 
reflect the impact of the momentum on the scale, the results of TS and SS terms are 
multiplied by a factor of 1 — e^^^^ which suppresses the low p contributions. 

We can predict Ds production if we substitute s for u in the calculations. And the 
RF of Ds is determined by Eq.® with a/b = mjmc ~ 3/10 (a = 3). The SPD of S'^{z) can 
be obtained with the FF of c quark splitting into K 22 1 (review the method in detail in 

n ri c !— .J V 

Refs.(5| and |6|). And then we can calculate the FFs of -Df ^ and Dg" which are necessary 
in the calculation of SS term. The predicted results of Dg production for 0-20% central 
collisions are shown in Figj2l The contributions of TS terms are higher than those of SS 
terms at pt > 0.4 GeV/c for or px > I GeV/c for Dg, because it is hard to produce 
c shower partons in the jets created in the heavy ion collisions due to the large mass of c 
quark. 
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FIG. 1: transverse momentum spectrum in central Au+Au collisions with the data at low 
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FIG. 2: Predicted Ds transverse momentum spectrum in central Au+Au collisions. 



IV. PRODUCTIONS OF 6-FLAVORED MESONS 



The most remarkable difference between the situation of 6-flavored mesons and charmed 
mesons is that all the FFs of b quark used in this paper are obtained at starting scale of 
Q = 2mb = 10 GeV/c. So we should try to get the SPDs initiated by b quark with the 



method shown in Refs.l5j and 



Q 



xD(x) 



dxi dx 



Xi Jo X2 



-{Sf{xi), Sf {x2)}RMixi,X2,x) 



where 



{S'!{x,),Sf{x2)} 



l[SKx,)sf{-^HSK- 



Xi 



-)Sfix2)] 



(17) 



(18) 



I- X2' 

reflecting the symmetrization of the leading parton momentum fraction. The RFs 
Rm{xi,X2,x) determined by Eq.(jl]) for 6-flavored mesons are denoted with a/b = md/m-^ ~ 
3/50 (a = 3) for B, a/b = nis/m-^ ~ 4.5/50 (a = 4.5) for Bg and a/b = rric/m-j^ ~ 3/10 



a = 3) for Be, respectively. With the FFs D{x) of b sphtting into vr, K [22'], D° |2l| and B 



23|, we can get the results of (/ = u{u) , d{d)) , S^, and S^. When we talk about FF 
of 6 —7- i?, we have in mind that there are four fragmentati on p rocesses b — ?■ B^, b ^ B^, 



b ^ B and b ^ B^ . Strong assumption are made in Ref. 



23| that the FFs of these four 



processes all coincide if the influence of the electroweak interactions is neglected. The results 
of these SPDs are shown in FigJS] 

Again, applying the SPDs 5*^ and 5*^ (S*^) to Eq. f|T7|l we can calculate the unknown FFs 
of b splitting into 6-fiavored meson Bg (Be) further. We use the Peterson fragmentation 



function 



24| 



D^(^) = ^T7T^^^IT-1^' (19) 
[(1 — x)^ + ex\^ 

which has been frequently applied in connection with the fragmentation of heavy quarks 
into their mesons. The two parameters are = 0.0201, e = 0.0177 for Bg, and = 0.0227, 
e = 0.0098 for Be, respectively. The reproduced FFs with the parametrized SPDs and the 
predicted FFs are shown in FigJH As the results have shown, all the fits are excellent. 

It's necessary to point out that b shower parton initiated by gluon is not included in this 
paper, since the FFs of gluon to 6-fiavored mesons at Q = 10 GeV/c are not known. So 
we calculate the shower parton distribution of b (Eq.Q) and the shower-shower {SS) term 
(Eq.( IT4l) ) with i = b. It's another difference in calculating 6-flavored mesons from that for 
charmed mesons. 

We present the transverse momentum spectra of 6-fiavored mesons B, Bg and Be for 
0-20% central collisions in Figj5]-Fig|71 The results show that the contributions of TT are 
much lower than TS and SS. The main reason is that the mass of b quark mi, = 5 GeV is 
much larger than those of light quarks and c quark, which leads the much smaller thermal 
parton distribution. To see this point clearly, we exhibit the thermal parton distribution 
densities T{p) of b, c with jb = 1c = 0.26 and other light quarks in Fig|8]for the purpose of 
comparisons. The distribution of b quark is almost eight orders magnitude lower than others. 
But this has not serious impact on the results of TS terms, because the contributions consist 
of two parts as shown in Eq. fll2p . Take B meson as an example, one part is Td{pi)Sj;{p2) 
and the other is Sd{pi)Ti{p2) ■ It's obvious that the major component of the TS term comes 
from the first part, which is not affected by the lower thermal parton distribution of b. So 
the TS contributions of B and Bs productions are higher than SS in the region of > 1 
GeV/c. While for Be, the combined effect of the lower thermal parton distributions of c and 
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b makes the 1~S contribution lower than SS. 

The fact that TS or SS contribution is much larger than TT contribution leads the 
different behavior of pt spectra for b-flavored mesons with that for charmed mesons. For 
D° and Ds the dominant contributions come from TT terms at low pt because the thermal 
parton distribution of c is higher than that of b. While for 6-fiavored mesons the primary 
behavior of the spectra are determined by TS or SS contribution. 
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FIG. 3: Shower parton distributions. 




V. MESON PRODUCTION FOR DIFFERENT CENTRALITIES 

In Ref.js], we have investigated the J/ip transverse momentum spectrum in central 
Au+Au collisions at midrapidity. In this paper we can predict the heavy flavored meson 
productions for non-central centralities. The fitted parameters involved can be obtained by 
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FIG. 5: Predicted B transverse momentum spectrum in central Au+Au collisions. 
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FIG. 6: Predicted Bg transverse momentum spectrum in central Au+Au collisions. 



fitting the experimental data of J ftp productions measured by PHENIX collaboration 25 |. 
Tlie flow velocity vt and heavy quark fugacity 7c in Eq.(l2]) (TT term) are the parameters 
adjusted to fit the experimental data at low px {pt < 2 GeV/c). We find that the data are 
fitted well with the changed 7c and the fixed vt = 0.3c as can be seen in FigjOl We suppose 




FIG. 7: Predicted Be transverse momentum spectrum in central Au+Au collisions. 
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FIG. 8: Thermal parton distributions. 

that fugacities for other quarks {u, d and s) are proportional to those for c quarks. 

When we calculate the TS and SS terms, /3L standing for the energy loss factor in the 
heavy ion collisions should be regarded as a function of centrality. And it is determined by 
the inclusive distribution of vr*^ for all centralities in Au+Au collisions pLfi]- The different 
values of 7 and /3L for different centralities are given in Table 1. For the hard parton 
distributions fi{k) is scaled by the number of binary collisions Ncoii- With the changed 
values of 7, /3L and fi{k) for the corresponding centralities we predict the spectra of heavy 
flavored mesons as shown in FigJTOl-FiglT^ for Dg., B, Bg and B^. for centrality bins 
0-20%, 20-40%, 40-60% and 60-92%. 

TABLE I: parameters for quark fugacities 7 and f3L for different centralities 





- 20% 20 


- 40% 40 


- 60 


% 60 - 92% 


7c = 76 


0.26 


0.19 


0.12 


0.06 


lu = Id 


1.00 


0.73 


0.46 


0.23 


Is 


0.80 


0.58 


0.37 


0.18 


(3L 


2.39 


1.44 


1.03 


0.78 



VI. CONCLUSIONS 

We complete the SPDs initiated by the hard partons of heavy quarks c and b by making 
use of the FFs in the framework of the recombination model. With the parametrized SPDs, 
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FIG. 9: J/ij) transverse momentum spectrum in Au+Au collisions for different centralities. The 



experimental data are from Ref. 




FIG. 10: Predicted transverse momentum spectrum in Au+Au collisions for different centrali- 
ties. 




FIG. 11: Predicted Dg transverse momentum spectrum in Au+Au collisions for different centrali- 
ties. 
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FIG. 12: Predicted B transverse momentum spectrum in Au+Au collisions for different centralities. 
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FIG. 13: Predicted Bs transverse momentum spectrum in Au+Au collisions for different centrali- 
ties. 

we reproduce the FFs of c or 6 splitting into charmed mesons or 6-flavored mesons and 
predict other unknown FFs. The meson spectrum in relativistic heavy ion collisions is the 
summation of the recombination of thermal-thermal, thermal-shower and shower-shower 
partons, where the shower-shower recombination equals to the subprocess of hard scattering 
followed by parton fragmentation. As the spectra of the heavy flavored mesons have shown, 
the contributions from thermal-shower recombination are higher than those from shower- 
shower recombination in the region of px > 1 GeV/c except for the results of Be- So the 
recombination of thermal- shower partons is important, of course, at every high pt, shower- 
shower recombination will be the dominant process for the heavy flavored meson production. 
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FIG. 14: Predicted Be transverse momentum spectrum in Au+Au collisions for different centrali- 
ties. 
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